猪基因共表达网络模块的构建及功能分析

doi:10.11843/j.issn.0366-6964.2017.11.023

摘要/Abstract

摘要：

旨在构建猪（Sus scrofa）的基因共表达网络模块，挖掘功能基因。基于权重基因共表达网络分析技术（Weighted gene co-expression network analysis，WGCNA），利用猪的多组织转录组测序数据，构建了24个猪的基因共表达模块，并对模块进行生物学过程、KEGG通路、疾病和组织特异性分析。结果显示，模块具有功能特异性，与器官的发育阶段、代谢通路有关，且同一模块中的信号通路间呈现调控关系；模块可富集与代谢、神经系统、癌症相关的疾病基因；模块可定位于特定组织。分析发现，钙调通路、PI3K-Akt通路、MAPK通路、泛素介导通路等在模块中呈现关键作用，催产素通路、昼夜节律通路内的相关基因对繁殖调控具有重要作用，SIX1、PRKAG3等基因的局部网络涉及肌肉发育。综上表明，在基因组水平上，构建和分析基因表达调控网络，可为探究中国地方猪品种特有的基因资源信息、丰富分子育种理论，提供新思路和新途径。

Abstract:

The present study aimed to construct gene co-expression network modules in pig, and explore functional genes. The 24 modules were constructed using RNA-seq dataset by Weighted Gene Co-expression Network Analysis, which were enriched in biological process, KEGG pathway, disease and tissue specificity. Modules showed significant functional specificity, and were correlated with specific development stages and metabolic pathways of organs. In addition, signaling pathways within the same module had regulation relationships with each other. Furthermore, modules enriched genes in diseases related to metabolic, nervous system and cancer, and also modules were located at the specific tissues. Moreover, calcium channel, PI3K-Akt signaling pathway, MAPK signaling pathway and ubiquitin mediating signaling pathway played key roles in different modules. The genes related to oxytocin and circadian entrainment pathways performed important regulatory functions in reproduction. SIX1 and PRKAG3 involved in development of muscles. Finally, the construction and analysis of gene co-expression regulation network can provide new ideas and methods to discover specific gene resource in Chinese local pig breeds and enrich the theories of molecular breeding.

中图分类号:

S828.2

魏凯, 张婷婷, 马磊. 猪基因共表达网络模块的构建及功能分析[J]. 畜牧兽医学报, 2017, 48(11): 2205-2215.

WEI Kai, ZHANG Ting-ting, MA Lei. Construction and Functional Analysis of Gene Co-expression Network Modules[J]. ACTA VETERINARIA ET ZOOTECHNICA SINICA, 2017, 48(11): 2205-2215.

参考文献

[1] TALUKDAR H A, ASL F H, JAIN R K, et al. Cross-tissue regulatory gene networks in coronary artery disease[J]. Cell Syst, 2016, 2(3):196-208.
[2] ALBERT R. Scale-free networks in cell biology[J]. J Cell Sci, 2005, 118(Pt 21):4947-4957.
[3] WEIRAUCH M T. Gene coexpression networks for the analysis of DNA microarray data[M]//DEHMER M, EMMERT-STREIB F, GRABER A, et al. Applied Statistics for Network Biology:Methods in Systems Biology. Weinhiem:Wiley-VCH Verlag GmbH & Co. KGaA, 2011:215-250.
[4] BARABÁSI A L, OLTVAI Z N. Network biology:Understanding the cell's functional organization[J]. Nat Rev Genet, 2004, 5(2):101-113.
[5] HARTWELL L H, HOPFIELD J J, LEIBLER S, et al. From molecular to modular cell biology[J]. Nature, 1999, 402(6761 Suppl):C47-C52.
[6] STUART J M, SEGAL E, KOLLER D, et al. A gene-coexpression network for global discovery of conserved genetic modules[J]. Science, 2003, 302(5643):249-255.
[7] KOGELMAN L J A, CIRERA S, ZHERNAKOVA D V, et al. Identification of co-expression gene networks, regulatory genes and pathways for obesity based on adipose tissue RNA Sequencing in a porcine model[J]. BMC Med Genomics, 2014, 7(1):57-73.
[8] JOSHI A, DE SMET R, MARCHAL K, et al. Module networks revisited:Computational assessment and prioritization of model predictions[J]. Bioinformatics, 2009, 25(4):490-496.
[9] KADARMIDEEN H N, WATSON-HAIGH N S. Building gene co-expression networks using transcriptomics data for systems biology investigations:Comparison of methods using microarray data[J]. Bioinformation, 2012, 8(18):855-861.
[10] ZHAO W, LANGFELDER P, FULLER T, et al. Weighted gene coexpression network analysis:State of the art[J]. J Biopharm Stat, 2010, 20(2):281-300.
[11] LI B S, TSOI L C, SWINDELL W R, et al. Transcriptome analysis of psoriasis in a large case-control sample:RNA-Seq provides insights into disease mechanisms[J]. J Invest Dermatol, 2014, 134(7):1828-1838.
[12] LEE J H, GAO C, PENG G D, et al. Analysis of transcriptome complexity through RNA sequencing in normal and failing murine hearts[J]. Circ Res, 2011, 109(12):1332-1341.
[13] DARLINGTON T M, EHRINGER M A, LARSON C, et al. Transcriptome analysis of inbred long sleep and inbred short sleep mice[J]. Genes Brain Behav, 2013, 12(2):263-274.
[14] FARAJZADEH L, HORNSH? J H, MOMENI J, et al. Pairwise comparisons of ten porcine tissues identify differential transcriptional regulation at the gene, isoform, promoter and transcription start site level[J]. Biochem Biophys Res Commun, 2013, 438(2):346-352.
[15] KODAMA Y, SHUMWAY M, LEINONEN R. The sequence read archive:Explosive growth of sequencing data[J]. Nucleic Acids Res, 2012, 40(D1):D54-D56.
[16] PATEL R K, JAIN M. NGS QC toolkit:A toolkit for quality control of next generation sequencing data[J]. PLoS One, 2012, 7(2):e30619.
[17] TRAPNELL C, PACHTER L, SALZBERG S L. TopHat:Discovering splice junctions with RNA-Seq[J]. Bioinformatics, 2009, 25(9):1105-1111.
[18] TRAPNELL C, WILLIAMS B A, PERTEA G. Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation[J]. Nat Biotechnol, 2010, 28(5):511-515.
[19] MORTAZAVI A, WILLIAMS B A, MCCUE K, et al. Mapping and quantifying mammalian transcriptomes by RNA-Seq[J]. Nat Methods, 2008, 5(7):621-628.
[20] LANGFELDER P, HORVATH S. WGCNA:An R package for weighted correlation network analysis[J]. BMC Bioinformatics, 2008, 9(1):559.
[21] HUANG D W, SHERMAN B T, LEMPICKI R A. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources[J]. Nat Protoc, 2009, 4(1):44-57.
[22] HUANG D W, SHERMAN B T, LEMPICKI R A. Bioinformatics enrichment tools:Paths toward the comprehensive functional analysis of large gene lists[J]. Nucleic Acids Res, 2009, 37(1):1-13.
[23] XIE C, MAO X Z, HUANG J J, et al. KOBAS 2.0:A web server for annotation and identification of enriched pathways and diseases[J]. Nucleic Acids Res, 2011, 39(Suppl. 2):W316-W322.
[24] SHANNON P, MARKIEL A, OZIER O, et al. Cytoscape:A software environment for integrated models of biomolecular interaction networks[J]. Genome Res, 2003, 13(11):2498-2504.
[25] NELSON D L, COX M M. Lehninger Principles of Biochemistry[M]. 5th ed. New York:W. H. Freeman and Company, 2008.
[26] BOILLY B, VERCOUTTER-EDOUART A S, HONDERMARCK H, et al. FGF signals for cell proliferation and migration through different pathways[J]. Cytokine Growth Factor Rev, 2000, 11(4):295-302.
[27] MATTSON M P. Neurodegenerative Disease[M]//MOOREN F C. Encyclopedia of Exercise Medicine in Health and Disease. Berlin Heidelberg:Springer, 2012.
[28] PIZZA V, AGRESTA A, D'ACUNTO C W, et al. Neuroinflamm-aging and neurodegenerative diseases:An overview[J]. CNS Neurol Disord Drug Targets, 2011, 10(5):621-634.
[29] KIRBY D M, SALEMI R, SUGIANA C, et al. NDUFS6 mutations are a novel cause of lethal neonatal mitochondrial complex I deficiency[J]. J Clin Invest, 2004, 114(6):837-845.
[30] BRANDT U. Energy converting NADH:Quinone oxidoreductase (complex I)[J]. Annu Rev Biochem, 2006, 75:69-92.
[31] D'HAESELEER P, LIANG S D, SOMOGYI R. Genetic network inference:From co-expression clustering to reverse engineering[J]. Bioinformatics, 2000, 16(8):707-726.
[32] DAS S, MEHER P K, RAI A, et al. Statistical approaches for gene selection, hub gene identification and module interaction in gene co-expression network analysis:An application to aluminum stress in soybean (Glycine max L.)[J]. PLoS One, 2017, 12(1):e0169605.
[33] WU W J, HUANG R H, WU Q H, et al. The role of Six1 in the genesis of muscle cell and skeletal muscle development[J]. Int J Biol Sci, 2014, 10(9):983-989.
[34] 陶新, 邓波, 门小明, 等. 不同猪种肉质相关基因Hal、RN和FTO的多态性研究[J]. 畜牧兽医学报, 2012, 43(5):676-683.
TAO X, DENG B, MEN X M, et al. Research on polymorphism of Hal, RN and FTO genes related to meat quality in different pig breeds[J]. Acta Veterinaria et Zootechnica Sinica, 2012, 43(5):676-683. (in Chinese)
[35] MILAN D, JEON J T, LOOFT C et al. A mutation in PRKAG3 associated with excess glycogen content in pig skeletal muscle[J]. Science, 2000, 288(5469):1248-1251.
[36] KALIA L V, LANG A E. Parkinson's disease[J]. Lancet, 2015, 386(9996):896-912.
[37] POEWE W. Psychosis in Parkinson's disease[J]. Mov Disord, 2003, 18(S6):80-87.
[38] GRINEVICH V, KNOBLOCH-BOLLMANN H S, ELIAVA M, et al. Assembling the puzzle:Pathways of oxytocin signaling in the brain[J]. Biol Psychiatry, 2015, 79(3):155-164.
[39] CHATTERJEE O, PATIL K, SAHU A, et al. An overview of the oxytocin-oxytocin receptor signaling network[J]. J Cell Commun Signal, 2016, 10(4):355-360.
[40] SAKAI T, ISHIDA N. Circadian rhythms of female mating activity governed by clock genes in Drosophila[J]. Proc Natl Acad Sci U S A, 2001, 98(16):9221-9225.